Topological insulator in the core of the superconducting vortex in graphene 
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The core of the vortex in a general superconducting order parameter in graphene is argued to be 
ordered, with the possible local order parameters forming the algebra 17(1) x CZ(3). A sufficiently 
strong Zeeman coupling of the magnetic field of the vortex to the electron spin breaks the degeneracy 
in the core in favor of the anomalous quantum Hall state. I consider a variety of superconducting 
condensates on the honeycomb lattice and demonstrate the surprising universality of this result. A 
way to experimentally determine the outcome of the possible competition between different types 
of orders in the core is proposed. 



Besides the usual metallic state Dirac quasiparti- 
cles in graphene may, at least in principle, form a mul- 
titude of insulating phases [1]. The insulators should be 
understood as different types of dynamically or externally 
generated masses of Dirac fermions, which preserve the 
quasi-relativistic invariance of the low-energy theory, but 
break some of the space-time symmetries of the Hamil- 
tonian 3. Yet another type of ordered phase of elec- 
trons on honeycomb lattice is the superconductor, which 
breaks the U(l) symmetry associated with the conserva- 
tion of the particle number The superconducting 
state should be inducible in graphene by the proximity 
effect for example [?! . At the level of the continuum Dirac 
Hamiltonian, superconducting states in graphene may be 
distinguished by their symmetry with respect to the in- 
version of the two inequivalent Dirac points, as the s and 
p waves, and with respect to the translational invariance, 
as uniform and not. The variety of such superconduct- 
ing states notwithstanding, I show here that the core of 
the superconducting vortex in undoped graphene is al- 
ways ordered, and that, at least in the simplest model, 
always in the same way: the core of the vortex is the 
Haldane-Kane-Melle (HKM) topological insulator 

The above result is obtained within two complemen- 
tary Dirac representations of the Bogoliubov-de Gennes 
(BdG) quasiparticle Hamiltonian, in which it arises in 
related, but distinct ways. The problem of the internal 
structure of the superconducting vortex in graphene may 
be mapped onto the general Dirac Hamiltonian in pres- 
ence of two-component mass-term with a twist. It is well 
known that the single zero-energy state in this situation 
would imply a quantum number fractionalization [lo| . 
In graphene, however, the spin- 1/2 of the electrons dic- 
tates two copies of Dirac fermions, and the concomitant 
doubling of the number of zero-energy states restores the 
usual quantization rule. The same zero-energy states, 
nevertheless, are responsible for another nontrivial prop- 
erty of the vortex core: it always harbors a finite local 
order parameter (OP) of some kind. An early example 
of this effect was provided in ref. [ll| , where the vortex 



configuration of the in-plane Neel OP on honeycomb lat- 
tice was shown to induce an out-of-plane component of 
the same OP. More recently, the superconducting vortex 
was argued to exhibit charge-density waves (CDW) and 
bond-density waves (BDW) in the core [l3 • Here it will 
be shown that in general the vortex of unit vorticity in the 
mass of the two-copy four-component Dirac Hamiltonian 
allows four order parameters, which close U{1) x CZ(3), 
where CZ(3) stands for the three-dimensional Clifford al- 
gebra. The members of the algebra depend on the type of 
the underlying order; nevertheless, in a superconductor 
of arbitrary symmetry, if the magnetic field in the vortex 
core is sufficiently strong the HKM anomalous quantum 
Hall state is universally the ground state. The effect is 
due to the Zeeman coupling of the magnetic field to the 
electron spin, which, quite unexpectedly, always selects 
the above state out of a wide variety of competing pos- 
sibilities. 

The HKM insulator is interesting in its own right, due 
to its quantized spin Hall effect, and the accompanying 
structure of the edge states [9|] . Such a topologically non- 
trivial ground state is favored by the spin-orbit interac- 
tions in graphene, which however are well known to be 
way too weak to lead to an observable effect. It was 
recently proposed that this ground state may also result 
from the electron-electron Coulomb repulsion in the pres- 
ence of a bulge in the graphene sheet [11] . The present 
considerations suggest that placing graphene on top of a 
type-II superconductor in a mixed state may provide an 
alternative, albeit local, realization of this elusive state 
of matter. 

The announced result will be first derived within 
the rotationally invariant representation of the paring 
Hamiltonian. Graphene exhibits gapless excitations 
near two inequivalent Dirac points at iif, with K =^ 
(1, l/-\/3)(27r/a-\/3), for example, where a is the lattice 
spacing 2\. One may form an eight- component Dirac- 
Nambu fcrmion as '^^ = (^+: ^L), where 
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"^lik.uj) = {ul{k,uj),vl{k,uj),(TU^„{-k,-uj),av^^{-k,-uj)), (1) 
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where cr = ± labels the projection of the electron spin 
along the z-axis, and k — K + p, with \p\ <^ K . 
and Va are the Grassmann variables corresponding to 
two triangular sublattices of the honeycomb lattice. The 
imaginary-time Lagrangian for the excitations near the 
Dirac points in this representation becomes 

L^-^\x,T){dr + Ho)-^{x,r), (2) 

where the single-particle Hamiltonian is i7o = I2 ® 
iloliVi, with 70 = 0-3(8) 0-3, 71 = l2'S)cr2, and 72 = /2<8)o-i, 
and i — 1,2. We may then choose 73 = ai (8) 173 and 
75 — (8 (73. {h, f?} is the Pauli basis of two-component 
matrices. 

The Lagrangian is invariant under the global U (4) 
symmetry generated by the 16 matrices {l2,<j} <E) 
{^4,73,75,735}, where 735 = 17375 = CT3 8) h, and I4 is 
the four-dimensional unit matrix. The SU{2) subgroup 
generated by the S — a (E) I4 is the group of rotations of 
the electron spin, and is, of course, exact. This is due 
to the factor "cr" in the right half of the Dirac-Nambu 
fermion in Eq. (1), which ensures that the rotation may 
be accomplished by the matrix multiplication from the 
left. The above representation is in this sense manifestly 
rotationally invariant. It will also prove helpful to recog- 
nize N = I2 (El 735 as the number operator. 

There are also 16 linearly independent Hermi- 
tian matrices that anticommute with Hq: {12,(7} (8 
{70, *7o73, *7o75, *7i72}- The addition of any of these 
matrices to the Hamiltonian Hq would gap the Dirac 
spectrum. The matrices that correspond to a super- 
conducting orders are those that do not commute with 
the particle- number. It is easy to see that (5'^(/2 81 
«7o73)*) = ^eA,, and (*t(/2 ^7075)*) = /mA,, 
where A^ represents the spatially uniform, complex s- 
wave superconducting OP. Consider then the BdG pair- 
ing Hamiltonian in the presence of the vortex in the un- 
derlying s-wave superconducting ground state: 

HBdG = Hq + I a, (r) 1/2 8) [no (73 sin6' + 75 cos 6*)], (3) 

where (r, 9) are the polar coordinates in the graphene 
plane. |As(r — >■ oo)| = const., and otherwise arbi- 
trary. The magnetic field in the vortex will be in- 
cluded shortly. An index theorem [3| guarantees that 
the spectrum of the block-diagonal Hamiltonian HbcIG 
contains two states with zero energy: ^'01 = ("00,0) 
and 2 — (O,'0o), where -00 is the rotationally invari- 
ant, four-component bound state [l3|- The Hilbert space 
Hq spanned by these two zero-energy states is invari- 
ant under all operators that commute or anticommute 
with HBdG [il], with four falling into the latter category: 
{I2, 17} 81 7o- These are at the center of our study. 

On the other hand, the expectation value of any trace- 
less single-particle operator M that anticommutes with 
the Dirac Hamiltonian is easily seen to derive entirely 
from the zero-energy states TT]. At T = 0: 

(*tM*) = - J2 )i'Ux)M^Jo.A^), (4) 

i,oc ijUnoc 



where {ipo^i{x)} is a basis in Hq. Since "Hq is two- 
dimensional, the four anticommuting operators from 
above in Hq reduce to the familiar Pauli matrices {I2, (?}. 
In full analogy to the spin- 1/2 problem [Tsj , any two or- 
thogonal states in Hq are then the -1-1 and —1 eigenstates 
of (n • (?) 8" 7o for some unit vector n , and the -1-1 eigen- 
states of I2 8) 7o. Equation (4) then implies that: a) 
when one state is occupied and the second state empty, 

. a 8) 70)*) = -0(5 -00, and (^'^(/z 8) 7o)*> = 0, 
and b) when both states are occupied or both empty 

(^'■|'(/2 8)70)^') = ±0(500, with (4''l'(n-CT8)7o)*) = for 
all h. 

It is straightforward to rewrite these fermion bilinears 
in terms of the original electrons. For example 

^''^(^2 1870)* = uI{x,t)u„{x,t) - v'l{x,T)Va{x,T). (5) 

The corresponding average represents therefore the 
CDW. Similarly, 

^'t(fc,a;)(cr3 <»7o)*(^,'^) = (7[ul{k,Lu)u^{k,uj) (6) 
-ul{-k, -ijj)ua{-k, -uj)] - [u ^ v], 

and may be recognized as the z-component of the vec- 
tor OP for the HKM topological insulator [1, Q. The 
ground state with a finite expectation value of this oper- 
ator would break the rotational invariance, the sublattice 
exchange symmetry, and, most importantly, the time re- 
versal for each spin projection separately. 

Let us take the effect of the magnetic field of the vortex 
into account next. Assume the field to be along the z- 
direction, orthogonal to the graphene plane. First, the 
sole orbital effect of the localized magnetic field is the 
change of the form of the zero-energy states [l^. This 
follows from the observation that in the Dirac-Nambu 
representation the magnetic field enters the Hamiltonian 
in Eq. (3) by modifying only the Hq term as 

Hq Hq [A] ^h® iJQJ^ {p^ - A^SS ) = (7) 
g-x(^)-f2(870^pg-x(5)-f2®70 

where the magnetic field is B = CijdiAj = d^x- (In 
our units h = e = c ^ 1.) Since the matrix I2 8" 70 
anticommutes with the mass term in Eq. (3), and acts 
as the unit operator within Hq, the zero-energy states 
with and without the magnetic field differ only in the 
factor exp[x(a;)]. More important turns out to be the 
Zeeman effect on the electron spin. With the Zeeman 
term the Hamiltonian becomes HBdG + Hz, where 

Hz = 9B{x){a3<El4), (8) 

and g ~ 2 for the electrons in graphene. One may con- 
sider Hz as a weak perturbation to the spectrum of HBdG 
[l7| . The commutation relations with the four operators 
that anticommute with HBdG imply that, within Hq, Hz 
is proportional to 173 8) 70, and therefore it splits the -1-1 
and —1 eigenstates of this particular OP. Equation (4) 
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then imphes that 

{'^^{as^jom ^\Exp[x{x)~ f A(r')dr']|2, (9) 

whereas the averages of the other three OPs vanish. As 
the total flux of the magnetic field is hc/2e, we have 
x{x) ~ (l/2)ln|x| at distances beyond the magnetic 
field's penetration depth, and the OP is exponentially 
localized within the coherence length ^ = 1/As(r = oo). 

Interestingly, the same conclusion follows even for the 
superconductor with a different symmetry. Consider the 
BdG Hamiltonian for a vortex in the p-wave state: 

HBdG = i?o + |Ap(r)|a3® [i7o(73sin6l + 75cos6')], (10) 

which breaks the spin-rotational symmetry. This super- 
conducting OP is odd under the exchange of the Dirac 
points, and it is energetically preferable for strong next- 
nearest neighbor attraction between electrons [5|. The 



T is the time-reversal operator, consisting of spin, mo- 
mentum, and frequency reversal, and particle-hole oper- 
ator exchange. The Lagrangian can still be written in 
the form in Eq. (2), but with the Hamiltonian Hq re- 
placed with Hq = [h ® i7o7i)pi + (fs ® «7o72)p2, and 
with 7o = /2'X'0'3, 7i = (^z®<^2, 72 = h®'^i^ 73 = o'i(g)(72, 
and 75 = 0-2 (8) (T2 [ill- 

Once the number operator in this representation 
is recognized as TV = as (g) J4, the set of six- 
teen matrices that anticommute with Hq may be di- 
vided into those corresponding to the superconducting 
({o-i,o'2} (8) {7i:«7o72,ni75,«7i73}) and the insulating 
OPs ({<T3,/2}(8){7o, j7o73,«7o75,«7i72})- In revealing the 
orders corresponding to these fermion bilinears it is useful 
to discern P = a-i® 735 as the generator of translations, 
and Is = I2® «7i75 as the K ^ —K inversion operator 
0. This immediately implies that the averages ($^Af<I>) 
represent: a) uniform s wave, for M = {a\^a2\ ® *7i75, 
b) nonuniform s wave, for M = {171,(72} ® 17072 , c) uni- 
form p wave, for M = {(7i,a'2} ® 47173, and d) nonuni- 
form p wave superconductor, for M — {cri, (72} ^ 7i- Ex- 
plicit rewriting of the above bilinears shows that the non- 
uniform superconducting states have the 2K periodicity 
of the OP, and thus represent the LOFF phases [l^ with 
the Kekule texture. Similarly, the insulating states are: 
a) CDW, for M = 0-3 (g) 70, b) spin-scalar Kekule BDW 
|20| . for M = {(T3 (g) 17073, 12 ® *7o75}, c) z-component of 
the HKM, for M = I2® 17172- Exchanging 0-3 and I2 in 
the above further yields, a) Neel, b) spin-vector BDW, 
and c) spin-scalar Haldane's OP, respectively. 

Let us reconsider the Hamiltonian with the vortex in 



U{1) X CZ(3) algebra of OPs in the core is now different: 
U{1) = {h X 70} and C/(3) = {0-3 g) 70, cti ® 47172, 0-2 ® 
47172}- The last two matrices correspond to the x and 
y components of the Neel OP. The orbital effect of the 
magnetic field in resolving the degeneracy of the zero- 
energy states is still null. Hz, on the other hand, within 
"Ho is again proportional to 173 (g 70, as dictated by the 
commutation relation of the Hz with the above set of 
OPs. This implies the same result in Eq. (9) follows for 
the p-wave vortex as well. 

Let us now construct a different representation of the 
problem in which the manifest rotational invariance of 
the Lagrangian in Eq. (2) is sacrificed [31 . Doing so 
will enable us to include the non-uniform insulating and 
superconducting states, and make contact with the recent 
work of Ghaemi et al. [12]. Consider a different eight- 
component Dirac fermion ~ (<i>^, (T((Ti (g) 12)$+)^), 
where the Dirac field is now rather standard 



(11) 
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the s-wave superconducting OP. In this representation it 
takes the form of 

HBdG = Ha + \Asir)\{ai sinO + cts cos6l) (g ^71 75. (12) 

The U{1) X CZ(3) OP algebra consists now of Cl{3) = 
{a'3g) 70, 0-3 '8)47073,^28)47075}, andC/(l) = {128)47172}- 
The CDW and BDW belong to the Cl{3) part, in agree- 
ment with the observation in ref. . In this representa- 
tion the Zeeman term becomes proportional to the unit 
matrix /2 gi £4. The Zeeman effect of the magnetic field 
is thus to shift both zero-energy states equally, so that 
both are above, or below, the Fermi level [21|. From Eq. 
(4) we see that it is the OP belonging to the U{1) factor 
that develops an expectation value in this case, which is 
precisely the HKM state again. 

Within the ^-representation it is evident that the 
above result holds irrespectively of the symmetry of the 
superconducting OP. For each one of the four possible 
choices the C/(l) part of the U{1) x CZ(3) algebra of the 
core states will always be the HKM, since the matrix 
47172 anticommutes with all the superconducting states, 
while at the same time it commutes with all the insulat- 
ing states, which are the exclusive members of the Cl{3). 
Since the Zeeman term is proportional to the unity oper- 
ator only the universal U{1) part of the algebra develops 
the expectation value for any symmetry of the OP. 

There is always the C/(l) x Cl{3) algebra of possi- 
ble OPs in the core, irrespectively of the type of or- 
der supporting the vortex. For example, in the 
representation a vortex in the Neel OP in the x-y plane, 
~ (ai C08 6 + a2 sin 0)g) 47172, in the core may exhibit the 
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<^>l{q,uj) = {ul{K + q,uj),vl{K + q,uj),ul{-K + q,uj),vl{-K + q,uj)). 
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p-wave superconductor, the CDW. and the z-component 
of the same Neel OP {a^ (g) «7i72) The most general 
Dirac Hamiltonian with a vortex can be written as 



Hgen = Mipi + M2P2 + O1M3 + O2M4 



(13) 



where (01,02) = |o(r)|(cos 0, sin0) is the OP, and {Mi} 
are four arbitrary eight-dimensional Hermitian matrices 
satisfying [Mi, Mj]+ = 2Sij. Since {Mi} form a represen- 
tation of the Clifford algebra CI (A), by a unitary trans- 
formation it can be reduced to a block-diagonal form 
Mi = tti ® ^i, where and f3i are four-dimensional 
[221 • Since {a^} and {Pi} then provide two (in general, 
different) Hermitian four-dimensional representations of 
C/(4), they are equivalent [l^. Any Hg^n may therefore 
be transformed into the block-diagonal form of Hsdo as 
given in Eq. (3), with the same four-dimensional matrix 
operator in both blocks. The only four operators that an- 
ticommute with such a Hamiltonian are obviously I2 ^ 70 
and fT (g) 7o, which form the algebra U{1) x Cl{3). 

If the superconducting OP in graphene is induced by a 
layer of a type-H superconductor laid on top, the Zee- 
man energy would roughly be « TJxkelvin, if the 
magnetic field H is expressed in tesla. The leading ef- 
fect of the graphene lattice, for example, comes from the 
^ 2K Fourier component of the vortex configuration, 
which in ^-representations is proportional to the matrix 
ai(E)l4. It is straightforward to check using our formalism 
that this perturbation favors the scalar BDW in the core, 
in agreement with the result of direct lattice calculation 
[1^ . The magnitude of this splitting should crudely be 
Cdw ~ \A\{a/Q w (102a/^)2kelvin, where a w lA is the 
lattice spacing. In the NbsGe with Tc ~ 23-?^ for ex- 
ample, edw ~ . The upper critical field, however, is 



~ 40 Tesla, so at the magnetic fields of few tesla the or- 
der in the core should predominantly be the HKM. It is 
also possible to increase the HKM component of the OP 
by adding a component to the magnetic field parallel to 
the superconducting layer so not to suppress the conden- 
sate. A distinct sign of the HKM OP in the core would 
then be the increase of the gap in the local density of 
states with the total magnetic field at a finite tempera- 
ture. This follows from the generalization of Eq. (4) to 
finite temperatures: 



(OP) = ^(tanh(i 



2kuT 



^'^^ ) + stanh(- 



2kBT 



■)), 



(14) 

with s = 1 for the HKM, and s = — 1 for a linear com- 
bination of the CDW and BDW. Assuming and edw 
to differ by a factor of 5 or more one OP is by an order 
of magnitude larger than the other over a wide range of 
temperatures. If the dominant order is any density wave, 
for e^u! 3> f-z 1 the OP would in contrast decrease with the 
magnetic field. 

It may also be useful to note that the competing or- 
ders behave differently under the changes of the signs of 
vorticity and/or magnetic field. For example, if the sign 
of both the vorticity and the magnetic field is reversed 
the BDW changes into its complementary pattern [l3| 
whereas the HKM state remains invariant. A probe sen- 
sitive to the sign of the gap would thus help distinguish 
between the different core states. 
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